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Abstract 

Motivated by the dark energy issue, the one-loop quantization approach for a class of 
relativistic higher order theories is discussed in some detail. A specific F(R, P, Q) gravity 

rvq ■ model at the one-loop level in a de Sitter universe is investigated, extending the similar 

program developed for the case of F(R) gravity. The stability conditions under arbitrary 

^D • perturbations are derived, 
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1 Introduction 

It is well known that recent astrophysical data indicate that our universe is currently in a phase 
of accelerated expansion. This is one of the most important achievement in cosmology. The 
origin of this observation is substantially not completely understood and the related issue is 
called the dark energy problem. 

Several possible explanations have been proposed in literature, among them one of the most 
popular is based on the use of gravitational modified models, the simplest one being Einstein 
gravity plus the inclusion of a small and positive cosmological constant, and this model works 
quite well but, having however some drawbacks (see, for example [HO [3] and reference therein). 

Roughly, the idea is that Einstein gravity is only an approximate low energy contribution, 
and additional terms depending on quadratic curvature invariants should be included. The 
idea is quite old, one of the first proposal was contained in [4J, where quantum R 2 gravity 
modifications were investigated (for a review, see [5]). The inclusion of general higher order 
contributions is also important for another aspect, since, sometimes they give extra terms which 
may also realize the early time inflation [6]. 

In previous papers [3 El [10] , f(R) gravity models and a non local Gauss-Bonnet gravity 
model at the one-loop level in a de Sitter background have been investigated. A similar pro- 
gram for the case of pure Einstein gravity was initiated in refs. |11[ [12| [13] (see also [141 I15j). 
Furthermore, such approach also suggests a possible way of investigating the cosmological con- 
stant issue [13]. Hence, the study of one-loop generalized modified gravity is a natural step to 
be undertaken for the completion of such a program, keeping always in mind, however, that a 
consistent quantum gravity theory is not available yet. 

Making use of generalized zeta- functions regularization (see, for instance [TU j IT7 1 \TE \ [TU | |2"U] ) . 
one may evaluate the one-loop effective action and then study the possibility of stabilization of 
the de Sitter background by quantum effects. Recall that in the one-loop approximation, the 
theory can be conveniently described by the (Euclidean) one- loop partition function (see [5]). 
For example, in the simplest case of a scalar field, one has 

Z = e-W ( D<j> e~S WW = e" r W . (1.1) 



Here I[(j> c ] is the classical action, evaluated on the background field cp c , while V is the one-loop 
effective action, which can be related to the determinant of the fluctuation operator L by 

T = -lnZ = / + -lndet — , (1.2) 

{j? being a renormalization parameter, which appears for dimensional reasons. Of course, in 
dealing with gauge theories, one needs a gauge braking term and the related F-P ghost contri- 
bution. 

The functional determinant may formally expressed by 

In det ^ = " l^ dt r 1 Tr e~ tL/ ^ . (1.3) 

M Jo 

Here the heat trace Tre - plays a preeminent role. In fact, for a second-order, elliptic non 
negative differential operator L in a boundaryless compact d-dimensional manifold, one has the 



small-t asymptotic heat trace expansion 

oo 

Tr e~ tL ~ Y. M L ) tJ ~ d/2 > ( L4 ) 

i=o 

where Aj(L) are the Seeley-deWitt coefficients J2TJ E2]- As a result, the expression (|1.3p is 
divergent and a regularization and renormalization are required. Zeta-function regularization 
may be implemented by [16] 



1 poo t e ~^ 1 

r(e) = I -2l dt WT7) T "" L " =I -Ts^ L ^' < 1J » 

where the zeta function associated with L is defined by 

1 r°° 
as\L) = — dtt^Tve-^, ((s\L/^)=^X(s\L). (1.6) 

r(s) Jo 

For a second order differential operator in 4-dimensions, the integral is convergent as soon as 
Res > 2. 

As a consequence, C( s l-^) is regular at the origin and one gets the well known result C(0|£) = 
Ai{V). This quantity is computable (see, for example, |23j). Furthermore, one may perform a 
Taylor expansion of the zeta function 

C( £ |L) = C(0|L) + C , (0) £ + O(e 2 ), (1.7) 

thus 

r(e)=/-> i ) + «^).o g ,H^f> +0 (e,. (1.8, 

As a result, one gets the one- loop divergences as well as finite contributions to the one-loop 
effective action in terms of the zeta function. With regard to this, a theory is one-loop renor- 
malizable as soon as the divergences can be cancelled in a consistent way by the renormalization 
of the bare coupling constants present in the classical action I. 

In this paper, we shall investigate modified generalized models, described by a Lagrangian 
density F(R,P,Q), where R is the Ricci scalar, and P = RijR 1 ^, and Q = Ruj TS R^ rs are 
quadratic curvature invariants. We do not include the Gauss-Bonnet topological invariant be- 
cause in four dimension, it can be expressed as G = R 2 — 4P + Q. 

After some considerations at classical level, the main part of the paper will deal with the 
one-loop evaluation of a particular but interesting F(R, P, Q) model on the de Sitter space, more 
exactly on its Euclidean version S{A). 

The paper ends with an application to the stability of the de Sitter space within the class of 
the modified gravitational models investigated. 

2 Linear perturbation of F(R, P, Q) model at classical level 

As warm up exercise, we shall begin with a some considerations at classical level. The equation 
of motion for general F(R, P, Q) model can be found in Ref. [23] and will not reported here. 
In fact, for our purposes, in this Section, it will be sufficient to consider only the trace of the 



equations of motion, which is trivial in Einstein gravity R = —k 2 T, but, for a general F(R, P, Q) 
model, reads 

A (3F' R + RF'p) + 2V i V j [{f' p + 2F' Q ) R ij ] - 2F + RF' R + 2 [PF' P + QF' Q ) = k 2 T . (2.1) 

Requiring R = Rq, constant and non negative, P = Pq, and Q = Qo constant, one has de Sitter 
existence condition in vacuum 



= 0. (2.2) 

R=R ,P=P ,Q=Q 



2F - RF' R - 2PF P - 2QF' Q 

As a particular but interesting model, let us make the choice 

F(R,P,Q) = f(R) + aP + bQ, (2.3) 

namely, a generic dependence on R, but only linear in the two quadratic invariants P and Q. 
With regard to this choice, as mentioned in the Introduction, the full quadratic case 

F(R,P,Q) = R-2A + cR 2 + aP + bQ. (2.4) 

namely a Einstein gravity with cosmological constant with the inclusion of curvature square 
terms is of particular interest and this model has been investigated in many papers, and was 
studied in the seminal paper [25jon the flat space. Note that in this particular case, we may take 
6 = 0, because the quadratic Gauss-Bonnet invariant G = R? — 4P + Q does not contribute to 
the equations of motion in four dimensions. Another interesting quadratic model is the Einstein 
plus conformal invariant quadratic term, i.e. 

F(R, P,Q) = R-2A + uj(^--2P + q\=R-2A + ujC ijrs C ijrs , (2.5) 

where Cij rs is the conformal invariant Weyl tensor. In the pure conformal case, one has the 
Weyl conformal gravity, and this is a quadratic model admitting exact black hole solutions (see 
for example [261EZ1E8]). 

Within the class of modified models (|2.3p . the dS existence condition becomes 



2/o-^o/o = 0. (2.6) 
and the trace equation in vacuum reads 

A (3/' + aR) + 2{a + 2b)V i V j R ij -2f + Rf' = 0. (2.7) 
Making use of contracted Bianchi Identity, 

ViVjR ij = ^AR, (2.8) 

one has 

A (3/' + 2{a + b)R) -2f + Rf' = 0. (2.9) 

Perturbing around dS space, namely R = Ro + 5R, one arrives at the perturbation equation 

- A5R + M$5R = 0, (2.10) 



in which the scalar degree of freedom effective mass reads 

M °-3^ + 2(a + 6)- (2 - U) 

Thus, Mq > is a necessary condition for the stability of the dS solution. In these F(R, P, Q) 
models, besides the massless graviton, there exists also a massive spin-two field, as we shall see 
in the next Section. 

In the particular case of a F(R,P,Q) = f(R) models, there is only the scalaron, and one 
recovers the well known condition for the dS stability (see, for example [291 GEH E] and references 
therein) 

f> 

'" >1. (2.12) 



Rofd 



3 Quantum field fluctuations around the maximally symmetric 
instantons 

In this Section we will discuss the one-loop quantization of the model on the a maximally 
symmetric space. Of course this should be considered only an effective approach (see, for 
instance [5]). To start with, we consider the Euclidean gravitational model described by the 
action 

I E [g} = - J d A x^F(R,P,Q) = - J d A x^[f(R) + aP + bQ], (3.1) 

with a, b are dimensionless (bare) parameters, the Newton constant G being included in the 
f(R) contribution. We assume the function F(R) to satisfy the condition (|2.2[) which ensures 
the existence of constant curvature solutions. This means that f(R) is not completely arbitrary, 
but it has to satisfy the equation 

/o -M = 0j (3.2) 

where here and in the following for the sake of simplicity we use the notation /o = f(Ro), 
f' G = f'(Ro) and so on. 

We are interested in the dS instanton S 4 with positive constant scalar curvature Rq. This is 
a maximally symmetric space having covariant conserved curvature tensors. Its metric may be 
written in the form 

ds 2 E = dr 2 (l - H 2 r 2 ) + _^ V) + r 2 dS 2 2 , (3.3) 

dS 2 being the metric of the two-dimensional sphere S 2 . The finite volume reads 

V (S*) = ^ , Ro = 12H 2 , (3.4) 

while Riemann and Ricci tensors are given by 

(o) _ Ro ( (o) (0) (0) (o)\ R (o) _ ^o (o) ,„ _v 
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Now let us consider small fluctuations around the maximally symmetric instanton. In the 
action (13,11) then we set 



9ij 



9ij • flij > 



9 



'j 



gV-h v +h**hi + 0(h d ), h = g^h ij , 



(3.6) 



,(°) 



where from now on gij = g\- is the metric of the maximally symmetric space and as usual, 
indices are lowered and raised by the means of such a metric. Up to second order in hij one has 



Va — > Va 



l + ~h+lh 2 -h ij h ij + 0{h 3 ) 

2 o 4 



(3.7) 



and 



R 



Ro 



Ro 



h + ViVih* 3 -Ah 



l v J 



+%h>% k 



I v t hV l h - i VthijVW + VihiVjh jk - i VjhikV'h^ , (3.8) 



where V& represents the covariant derivative in the unperturbed metric gij. More complicated 
expressions are obtained for the other invariants P, Q, but for our aim it is not necessary to 
write them explicitly. 

By performing a Taylor expansion of the Lagrangian around de Sitter metric, up to second 
order in hij, we get 



lE[g] ~ - / d 4 Xy/g 



hX 
F(R ,Po,Qo) + — + C 2 



(3.9) 



where £2 represents the second-order contribution and X = /q — Rq/q/2 vanishes when f(R) 
satisfies the de Sitter existence solution (J3.2J) . 

It is convenient to carry out the standard expansion of the tensor field hij in irreducible 
components PS], namely 



hj = h^ + V»£j + Vjii + ViVjO + - gijih - A a) 



(3.10) 



where a is the scalar component, while ^j and hij are the vector and tensor components with 
the properties 



ViC = 



Vih v = 



K = o 



(3.11) 



In terms of the irreducible components of the hij field, the Lagrangian density, disregarding 
total derivatives, becomes 



£2 



£hh + 2 Cher + C-aa + Cv +£-T: 



(3.12) 



where Chh, C-ha-, £-o<r represent the scalar contribution (a 2x 2 matrix), while Ly and Ct represent 
the vector and tensor contributions respectively. One has 



£ 



hi, 



h 



9fS A 2 3/^A aR A bRoA 3f{>R A , f>>R 2 f^Ro 



32 



32 



16 



X 3aA 2 36A 2 
+ 16 + 16 + 16 



16 
h, 



16 



:\2 



32 



(3.13) 



'ha 



9/^A 3 



+ 



3/^A 



32 ' 32 10 

3aA 3 36A 3 



^/^ A 2 - i-fHRlA + 



16 



16 



32" 

— aR A 2 bR A 2 

16 16 



32 



a, 



(3.14) 



(7 



//A 4 



9^A 
32 



_3^A3 
32 

3XA 2 
16 



+ ^/^oA 3 + ^«A 2 



32 



foRoA 2 



R XA 3aA 4 36A 4 



16 16 



1 
16' 



+ — — + -^aR A 3 + —bRnA 3 
16 



1 
16' 



(3.15) 



C 



V 



-XR + -XA 

8 2 



£,k , 



(3.16) 



£r 



h 1 -' 



+ 



f A f^Ro XaA 2 , ,. 2 , ai? A 6fl A aR 2 , 6ft 2 



24 



4 4 



+ 6A^ + 



24 



72 



+ 



36 



ha • (3.17) 



where A = g^VjVj is the Laplace-Beltrami operator in the unperturbed metric gij, which is 
a solution of field equations, but only if X = 0. We have written the above expansions around 
a maximally symmetric space which in principle could not be a solution. This means that the 
function f(R) can be arbitrary. 

As it is well known, invariance under diffeormorphisms renders the operator in the (h, a) 
sector not invertible. One needs a gauge fixing term and a corresponding ghost compensating 
term. Here we choose the harmonic gauge, that is 



Xj = -Vih t j -^V j h = 0, 

and the gauge fixing term 

C 9f = 2X l G ij x', G ij =~/g ij . 

The corresponding ghost Lagrangian reads [5] 
5 X k 



C 



yii 



B l Gi 



5 si 



-C J 



(3.18) 



(3.19) 



(3.20) 



where Ck and B^ are the ghost and anti-ghost vector fields respectively, while 5 x is the variation 
of the gauge condition due to an infinitesimal gauge transformation of the field. In this case it 
reads 



5 hi 



Vi£j + Vj£i 



bei 



gij A + Rij . 



(3.21) 



Neglecting total derivatives, one has 



L gh = B* 1 A + 



Re 



Cu . 



In irreducible components one finally obtains 



C 



of 



1 
2 



'' "(Al+f) ^+ 3 {hU + f)A lt rr 





P 2 r a , 9 /a R oV a 




£gh = 7 


[^(A 1 + f)c t + ^S(A --^)A„ £ ] 


i ghost irreducible components are defined by 


C k = C k + V k c, v k c k = o, 


B k = B k + V k b, V k B k = 0. 





(3.22) 



(3.23) 
(3.24) 

(3.25) 



4 One-loop effective action 

In order to compute the one-loop contributions to the effective action one has to consider the 
path integral for the bilinear part C = £2 + £> g f + £- g h of the total Lagrangian and take into 
account the Jacobian due to the change of variables with respect to the original ones. In this 
way one gets [HI IS] 



Z 



(i) 



-1/2 



(det Gy)- i/z J D[hij\D[C k ]D[B k ] exp (- j d 4 x^C 

(det Gij)~ 1/2 det J^ 1 det J 2 1/2 

DlhjDlh^Dl^DlajDlC^DlB^DlclDlb] exp 



d Xy/gC 



(4.1) 



where Ji and J2 are the Jacobians due to the change of variables in the ghost and tensor sectors 
respectively [13]. They read 



■h 



■h 



Ai + 



Ao + f )A„ 



(4.2) 



and the determinant of the operator Gy in this case is trivial. Due to the presence of curvature, 
the Euclidean gravitational action is not bounded from below, because arbitrary negative con- 
tributions can be induced on R, by conformal rescaling of the metric. For this reason, we have 
also used the Hawking prescription of integrating over imaginary scalar fields. Furthermore, the 
problem of the presence of additional zero modes introduced by the decomposition (|3.10p can 
be treated making use of the method presented in Ref . [13] . 

Now, a straightforward computation leads to the following off-shell one-loop contribution to 
the "partition function" 

-1/2 



-r« 



Z (D 



det 



A, + ^ + ^)(A„ + 5> + £]l„-l n -lJl; 



x det 



A + ^ 



Ro 

Al + T 



(4.3) 



where 



U 



2/o 



5Rofs 



2R (a + b) 



4(3# + 2(a + 6)] 



± 



[Wo ~ SRofo - 2^0(0 + b)} 2 - 8[3/£ + 2(o + b)][2X - RoW - RofH)] 



4(3/£ + 2(a + 6)] 



L 



A 2 + 



6f^ + Ro(a 



{2f + aR ) 2 + 16X(a + 46) 



12(a + 46) ^ 4(a + 46) 

A 0, A 1, A 2 being respectively the Laplacian operators acting on scalars, transverse vectors and 
transverse, traceless tensors and of course a + 46 7^ 0. 

The partition function in (|4.3p explicitly depends on the gauge parameter 7, but it is known 
that when one goes "on-shell", that is when one imposes the background metric gy to be a 
solution of the field equation, the one-loop partition function becomes gauge independent. In 
our case we have simply to perform the limit X — > obtaining 



7(1) 

on— shell 



det 



-A, 



x <^ det 



jo - Rofo 
3/q 7 + 2(a + 6) 

Rq 

6 



-1/2 



det 



-Ai-$ 



-A 2 + 



-A 2 + 



(26 - a)i? - 3/^ 



3(a + 46) 
As a consequence, the on-shell one-loop effective action reads 

r on-shell = Ie(9) + ^on-shell ' 



1/2 



-1/2 



(4.4) 



,(1) 

on— shell 



1 



log det 



-A + 



1 



log det 



1 



I* 1 



— log det 
2 B 



1 



M 



/o ~ go/p 
3Z6 7 + 2(a + 6) 

a -Ro 

-A 2 + — - 
6 



+- log det 



/' z 



-A 2 + 



(26 -a)R - 3/6" 
3(o + 46) 



(4.5) 



As usual an arbitrary renormalization parameter l/^ 2 has been introduced for dimensional 
reasons. 

The one-loop contribution r' 1 ^ to the effective action can be computed by by using zeta- 
function techniques. The eigenvalues of Laplacian operators on S 4 are explicitly known and so 
the determinant of all operators appearing in (|4,5p in principle can be calculated. We refer the 
interested reader to Ref. [7J, where all details of computation can be found. 



5 Discussion and conclusions 



We conclude the paper with several remarks. First, the one-loop effective action result is in 
agreement with a similar one given in [9] and in the limit a — > 0, 6 — > becomes identical to 
those in [7J, where only f(R) modified gravity has been considered. 



Both equations in (|4.3|) and (|4.4|) have been derived by assuming a + 46 ^ 0. As a useful 
check, we note that when a = —46, the original classical Lagrangian density can be written in 
the form 



f(R) + a P + bQ = f(R) + bG . 



f(R) = f(R) -bR 2 , G = R 2 -4P + Q. 



where G is the Gauss-Bonnet topological invariant which does not contribute to the field equa- 
tions. In such a case, as expected, the one-loop contribution to the partition function becomes 



Z 



(i) 

on— shell 



det 



-An + 




A RQ 

-A 2 + — - 
6 



-1/2 



which again is the result obtained for a pure f{R) modified gravity [7]. 

Another interesting particular case is the one in which the Einstein-Hilbert Lagrangian den- 
sity is modified by a term proportional to the the square of the Weyl tensor. Such a classical 
model has a de Sitter solution and in the absence of Einstein-Hilbert term, the on-shell, one-loop 
contribution to the partition function trivially vanishes, but if a = —46, 6 = u and, as in (|2.5p 



f(R) + aP + bQ = f(R) + ul^--2P + Q 



f(R) 



R 



16ttGn 



M 2 P R, 



Gn,Mp being respctively the Newton constant and Planck mass, the on-shell one- loop contri- 
bution reads 

-1/2 



7(1) 
^Weyl 



det -A 



Rq 

4 



1/2 



det 



A 2 + 



Ro_ 
6 



-A 2 + 



Rq 



In contrast with previous cases, here the contribution due to scalar components is vanishing. 

As an important application we discuss the stability of de Sitter space. To this aim we have 
to recall that the eigenvalues A n of Laplacian-Beltrami —A operators on S 4 have the form 



Xn ~ V2 



(n + vf 



a 



9r, 



ci(n + v) + c 3 (n + u) 3 , 



n 



0,1,2,. 



g n being the corresponding degeneracy and v, c\, C3 being dimensionless quantities, which depend 
on the operator one is dealing with. In particular one has 



-Ao 
-A 1 
-A 2 



UL — 45 


Cl = 
Cl = 


1 

12 ' 
9 
4 ' 


" 4 5 


Cl = 


125 
12 



C'3 
C3 
C3 



(5.1) 



We see that only the scalar Laplacian — Ao has a null eigenvalue, while the minimum eigenvalue 
of —A 1 is -Ro/4 and the minimum eigenvalue of — A 2 is 2Rq/3. 

In Eq. (|4.5p we are dealing with operators of the kind L = —A + M 2 and so, in order to 
have stability of de Sitter solution, we have to assume all eigenvalues of L to be positive. In this 
way we obtain restrictions on the function F(R, P, Q). 

Looking at (|4.5p we see that, independently on the classical action, a zero-mode is present 
coming from the Laplacian-like operator —A 1 — Rq/4. In principle, other zero modes may be 
present and all of these can be treated according to Ref. |13| . 



The other operators in (|4.5|) which could have vanishing or negative eigenvalues are 



2 Tl/tI 



L = -A + M£, M£ 



/o ~ Rofo 
3fg + 2(a + 6) ' 

L 2 = -A 2 + M> M 2 = Ro(2b-a)-3K 



3(a + 46) 



but in the case in which 



Ml > =► JL^/L > . 



1/2 , 2d ^ n >> (0+106)^0-3^ „ 



".Vf 2 ** 6 ^, " ' (5.2) 

t 2 -r 3^0 ^u = 3(a+4fe) ° " " 

In the particular cases in which Mq = and/or M% = —2/3Rq there are other zero-modes which 
have to be treated as the previous ones. 

In the interesting case f(R) = Mp(R—2A), linear in the curvature, the dS stability conditions 
become 



M 2 - M p { a + b > ° 

M " 2(a+b) » I Q + 46 > 

2 _ 4A(q+10b)-3M^ =^ S 2 

M 2 ~ 3 (a+4fe) > U . a + 106 > ^ 



a + 6 > 

a + 46<0 . (5.3) 

a + 106 < ^ 



4A 

We see that depending on the arbitrary parameters a, 6 the solution can be stable or unstable 
and in order to have a stable solution at least one of the two parameters has to be positive. In 
particular, in the special cases a = or 6 = one gets the stability conditions 

3Ml 
a> — £ > 0, 6 = 0, 

4A 

, 3Ml 

b> id >0 - °=°- 

In summary, here we have evaluated the one-loop effective action for a specific modified 
gravity model in de Sitter space. Generalized zeta regularization could be used to obtain a finite 
answer for the functional determinants in the effective action, what has proven to be a very 
convenient procedure. 

The important lesson to be drawn from this calculation, generalizing the previous program 
for one- loop Einstein gravity and f{R) modified gravity in the de Sitter background is that 
quantum corrections may to destabilize the classical de Sitter universe, as we have explicitly 
verified in the examples. 
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